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GEOMETRIC THEORY OF LATTICE VIBRATIONS 
AND SPECIFIC HEAT 

MIKHAIL SHUBIN AND TOSHIKAZU SUNADA 


Abstract. We discuss, from a geometric standpoint, the specific 
heat of a solid. This is a classical subject in solid state physics 
which dates back to a pioneering work by Einstein (1907) and its 
refinement by Debye (1912). Using a special quantization of crys¬ 
tal lattices and calculating the asymptotic of the integrated den¬ 
sity of states at the bottom of the spectrum, we obtain a rigorous 
derivation of the classical Debye law on the specific heat at low 
temperatures. The idea and method are taken from discrete geo¬ 
metric analysis which has been recently developed for the spectral 
geometry of crystal lattices. 


1. Introduction 

The primary purpose of this note is to discuss some dynamical prop¬ 
erties of solids, more specifically lattice vibrations in crystalline solids, 
from a geometric view point. In particular, we are concerned with a 
mathematically sound computation of the specific heat, a typical ther¬ 
modynamic quantity in solid state physics. The main idea is to employ 
a technique in discrete geometric analysis developed originally for the 
study of random walks on crystal lattices (uni, nn), and mathematical 
apparatus such as von Neumann trace and direct integrals which make 
the discussion more transparent than the existing ones. 

Theoretical computation of the specific heat at low temperature had 
been one of the central themes in quantum physics at the beginning 
of the last century (see jlH] and |lj for the history). The crucial point 
in the computation is to regard a solid as a crystal lattice realized 
periodically in the space R^. 

To explain what crystal lattices mean, let V be the set of constituent 
atoms in a solid, and : V —> R^ be the injective map representing an 
arrangement of atoms in equilibrium positions. Considering elements 
in V to be vertices, we join two elements in V by an (abstract) edge if 
they (as atoms) are bound by atomic forces, and extend to the set 
of edges as a piecewise linear map. We thus have a graph X = (V, E) 
realized in R^ where E denotes the set of all oriented edges. When 
<P(A) is invariant under the action of a lattice L C R^ by translations, 
the graph X (or its realization <P{X)) is said to be a crystal lattice. 

The inter-atomic forces allow the vibrations of atoms which involve 
small excursions from the equilibrium positions. We may describe the 
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vibration by the (linearized) equation of motion 


(fi 

dt'^ 


= Df, 


where D is a certain linear difference operator of the “second order” 
on X involving masses of atoms and inter-atomic forces in the coeffi¬ 
cients, and f = f(t, x) stands for the displacement from the equilibrium 
position (thus the position of an atom x at time t is ^{x) -|- f(t, x)). 

Following Planck’s idea on “energy quanta” originally applied to 
black-body radiation (1900), or according to the quantum mechanics 
founded by Heisenberg and Schrodinger (1925-26), physicists usually 
go forward as follows. 

(1) Consider the lattice vibration as an inhnite-dimensional system 
of harmonic oscillators. 

(2) Decompose the system into independent simple harmonic oscil¬ 
lators, and calculate the distribution of vibration frequencies. 

(3) Apply statistical mechanics to determine the macroscopic equi¬ 
librium state (the Gibbs state) of the quantized lattice vibration, and 
compute the internal energy U = U{T) (per unit cell) where T is the 
absolute temperature. Then the specific heat is given by 


C(T) 


dU 

&r' 


(See any textbook of solid state physics, for instance, [7j, P, or a review 
paper [2j, for the detail of this procedure in which one may see a daring 
manner of physicists to bring us effectively to the correct result. An 
exception is the book P by Born and Huang, which, in the second half, 
is written in a strictly deductive style.) 


This set routine leads us to the expression 


( 1 ) 


where 


U(T) = const + 



n^/x 

gh^IKT _ 


d(p(A), 


h = 2Tih = Planck constant, 

K = Boltzmann constant, 

and </5(A) is the (integrated) density of states satisfying the normaliza¬ 
tion condition 

d(p{\) = 3n [n = the number of atoms in a unit cell). 

The function (p(A) is closely related to the distribution of vibration 
frequencies (actually, \/X/27r represents the frequency parameter, and 
hence h\/X stands for “energy quanta”). As is easily seen, the behavior 
of the specific heat at low temperature relies heavily on the asymptotic 
behavior of (p(A) around A = 0. In the early stages of quantum physics. 
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however, physicists had no rigorous methods, in marked contrast to the 
case of the black-body radiation, to acquire precise information on (p(A) 
through the microscopic structure of a solid so that they were forced 
to make daring hypothesis on the shape of 9 ?(A) (indeed, it was in 1912 
when the discrete structure of solids were conhrmed by means of the 
diffraction of X-rays). 

The hrst substantial result was established by A. Einstein in 1907 
(two years after the publication of his three famous papers; [HI). He 
adopted a function </9(A) dehned by 


V7(A) 


0 (A < Ao) 
3n (A > Ao) 


with one characteristic frequency A/Ao/27r, or equivalently 

^ = 3n6{X - Ao) 

(i.e. replacing the vibration spectrum by a set of oscillators at a single 
frequency) to claim 

3„K(^y 

V KT J - l)^ 

In spite of the seemingly unrealistic model, this formula explains well 
not only the law of Dulong and Petit at high temperature : 



C{T) = 3nK, 


which had been known since 1791 in the setting of classical mechanics, 
but also the qualitative fact that, if T goes to zero, then so does C(T) 
as experimental results show (H. Nernst; 1910). Note, however, that 
Einstein’s formula for C{T) has exponential decay as T J, 0, which 
turns out to be quite incorrect. 

In 1912 (|S]), Debye proposed, without any knowledge of the lattice 
structure of crystals as in Einstein’s case, to take the function 


( 2 ) 


V7(A) 


f 0 (A < 0) 

< cqA^/^ (0 < a < \d) , 

[cqA^^ (A > Ad) 


with a suitable positive constant cq. The quantity Ad is taken so as to 
satisfy 

dip{\) = 3n. 

Debye’s distribution ip may be inferred from the intuitive observation 
that the continuum limit of a crystal lattice is a (uniform) elastic body, 
and that, in the limit, the density of states for lattice vibrations may 
be replaced by the one for elastic waves in a region of low frequencies. 
Actually, the constant Cq is determined in such a way that the function 
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is the integrated density states for elastic waves (see Section IT0|l . 
This view is natural because to one not aware of the atomic constitution 
of solids, a solid appears as an elastic continuum. In particular, it is 
deduced that, if the elastic body is isotropic, then the constant cq is 
given by 

6?!^ \cf cf / 

where 

Q = the longitudinal phase velocity, 

Ct = the transverse phase velocity, 

V = the volume of the unit cell 

(as a matter of fact, the crystalline solids are never isotropic, so that the 
constant Cq above should be replaced by an average of phase velocities 
over all directions of propagation; see jl] or Section ED- 

With his choice of the distribution (p, Debye gave the following neat 
formula for the specihc heat 

/ T \ 3 r^o/T 4 X 

C(T)=9nK(—) 2 

The quantity 0 _d is what we call the Debye temperature. It is dehned 
by 



In particular, it is found that, as T goes to zero, 

CD ^ 

2 /I 2 \ 

~ — H—(in the isotropic case). 

10 ' 

This is the T^-law which agrees well with the experimental data. 

Debye’s model m is still crude since the real shape of the distribution 
(p(A) turns out to be quite different from Q in the region of high 
frequencies. As far as the T^-law is concerned, however, we only need 
the asymptotic property 

(3) (^(A) ~ coA3/2 (A i 0). 

To establish this asymptotic behavior in a rigorous way, we shall ob¬ 
serve, on the one hand, that V3(A) is expressed as the von Neumann 
trace of the projection E{X) appearing in the spectral resolution of the 
operator D. On the other hand, we see that D is decomposed into a 
direct integral over the unitary character group of the lattice L (this 
corresponds to the decomposition of the lattice vibration into indepen¬ 
dent simple harmonic oscillators). In these discussions, a crucial fact is 
that D commutes with the action of the lattice on the crystal lattice. 
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A discrete analogue of trace formulae elucidates a close relation be¬ 
tween the function ip and the family of perturbed operators appearing 
in this direct integral (actually the formula is derived from this ob¬ 
servation). Up to this point, we do not require a detailed form of the 
operator D. To proceed further, we must impose a special condition 
on the matrix of atomic force constants which seems appropriate from 
both the nature of dynamics and the geometric view. This condition 
combined with a standard perturbation technique allows us to establish 
m without resorting to Debye’s continuum theory. 

Overall, we shall follow, in a slightly different fashion, the journey 
which physicists usually make to perform computation of the specihc 
heat. Besides some mathematical tools and ideas, the main difference 
is in the use of the terminology in graph theory which has a great 
advantage: not only allows it to avoid the redundancy of suffixes in 
the formulas, but also naturally brings geometric ideas at our disposal. 
This is the reason why we describe the crystal lattice as a graph. 

2. Crystal lattices 

We shall use the following notations. Given an oriented edge e E E 
of a graph X = (U, E), we put 

oe = the origin of e, 
te = the terminus of e, 
e = the inverse edge of e, 

and write, for x E V, 

Ex = {e E E] oe = x}. 

Let X = (U, E) be a crystal lattice realized in R^. We will always 
assume that the graph X is connected. (This assumption is natural 
if we wish X to imitate a single piece of a real crystal.) The lattice 
L acts freely on X by graph-automorphisms through the map <P. We 
write the action of L on U and E as (a, x) i—*• ax and (a, e) h-^ ae, 
respectively. Hence <P{ax) = ^(t) -|- a. We denote by Vq (resp. Eq) 
the quotient set of V (resp. E) by the action of L. Then a graph 
structure is induced on Xq = (Vq, Eq) in such a way that the crystal 
lattice X is an inhnite-fold abelian covering graph over Xq with the 
covering transformation group L. The set Vq is obviously finite. We 
assume that Eq is also hnite (thus we are treating the case of hnite 
range interaction). For the later purpose, we take a fundamental set T 
in V for the L-action. 

Assume that two atoms x and y are of the same kind if they are 
in the same L-orbit, that is, if there exists a E L such that y = ax. 
Denote by m{x) the mass of the atom x. As a function on V, m is 
L-invariant, and hence is regarded as a function on Vq. 
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We fix a unit cell, that is, a fundamental domain for the action of L 
on by translations (for instance, take a fundamental parallelotope). 
The number of atoms in a unit cell coincides with the number of vertices 
in Xq, which we denote by n. We write m(Vo) = which is 

the total mass of atoms in the unit cell. 

We denote by a ■ b the standard inner product on R^. Let L* be the 
dual lattice of L dehned by 

L* = {?7 G R^; ri-a eZ {a e L)}. 

The lattice 2nL* is what physicists usually call the reciprocal lattice. 
Thus 2717] {rj G R^) play the role of wavenumber vectors in the physical 
context. 

Examples (1) The cubic lattice X = (V, E) is a crystal lattice with 
V = Z^ C R^. Two vertices (mi,m 2 ,m 3 ), {ni,n2,n^) are joined by 
an edge if and only if ~ ''^*1 = 1- For the standard lattice 

L = Z^ acting on X in a natural manner, the quotient graph Xq is the 
3-bouquet graph consisting of a unique vertex with three loop edges. 

(2) The diamond lattice is defined as follows. Let 61 , 62,63 be the 
standard basis of R^, and let L be the lattice generated by 61 -I- 62 , 62 -|- 
63 , 63 - 1 - 61 . Then put 

l/ = Lu(L + (l/2,l/2,l/2)). 

Two vertices joined by an edge should have the forms 

(m 2 -h m 3 , m 3 mi, mi -h m 2 ), {n 2 + n 3 + ^,n 3 + ni + ni + n 2 + ^) 

satisfying one of the following conditions 

(i) mi = ni, m 2 = n2, m 3 = n^, 

(ii) mi = ni + 1 , m 2 = n 2 , m 3 = 

(hi) mi = ni, m 2 = n 2 + 1 , m 3 = 

(iv) mi = ni, m 2 =n 2 , m 3 = ns- 1 - 1 . 

It is easily checked that the quotient graph Xq by the action of L on 
X is the graph with two vertices joined by 4 multiple edges. 

3. The equation of motion 

The purpose of this section is to seek the form of the operator D. 
To this end, we impose several conditions on the dynamics of lattice 
vibrations. The discussion here is rather formal in the sense that we 
do not specify the domains of linear operators which we introduce. 

The lattice vibrations are supposed to be governed by a potential 
energy M<p(f) = u{<d> + i). Expand it for small / as 

n^(f) = n^(0) + G,p-i + ■ f H-, 
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where f ■ g = Yhx&v ' s(^)) cind K,p is a linear operator acting on 
displacements and satisfying 


( 4 ) 


■ g = f ■ K^g. 


The hrst order term ■ f must be zero, and ■ f > 0 because of the 
stability of the equilibrium positions. By omitting higher order terms 
as usual (i.e. neglecting a coupling produced by “anharmonic” terms), 
we write 

M<p(f) = M<f(0) + ■ f, 

and get the equation of motion 


( 5 ) 


m{x) — {x) = -grad u^, = -{K 4 ,t){x). 


Since the atomic force acts between two atoms x and y if and only if 
they are joined by an edge, we may write 


— {K,pi)(x) = A,p{e)i{te) + B,p{x)f{x), 


where A^{e) and B^{x) are linear transformations of R^. 

We impose a somewhat strong assumption on K^. We assume that, 
if is another realization obtained by a rigid motion of then the 
equation 

(72 f 

is equivalent to © provided that ^ + f = ^' + f'. This property leads 
to 

We write K = A = A,p, B = B^. 

Applying — ^') = 0 to the case + v, (n G R^), we get 

J2Me) + B{x) = 0, 

cGEx 

or equivalently 


-{Ki){x) = A(e)(f(te) - f(oe)), 

eGEx 

so that if we put 

^ ^ 41(e) (f(te) - f(oe)), 

^ ^ eeE^ 


then the equation of motion is given by 
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Applying also — ^') = 0 to the case U € SO {3), we 

have 

Y1 Me)Uv{e) = Me)v{e), 

eSEx gGEx 

where v(e) = ^(te) — ${oe) which may be regarded as a function on 
Eq due to the periodicity. Since {U — I; U E 50(3)} spans M 3 (R), 
the space of all real 3x3 matrices (see Lemma o for a proof), we 
conclude 

A(e)Tv{e) = 0 

e^Ex 

for every T G M 3 (R), or equivalently 


( 6 ) 


A{e) ® v{e) = 0 (as a tensor), 

CGEx 


that is 

^ A{e)ijv{e)k = 0 

CGEx 

for every 

We call A{e) the matrix of atomic force constants. From the nature 
of crystals, it is natural to assume that A{ae) = A{e) (a G L) (hence 
A{e) is regarded as a matrix-valued function on Eq). The symmetry 
condition ■ g = f ■ K,pg is equivalent to A(e) = *A(e). If A{e) is 
symmetric, i.e. A{e) = A{e), then the condition ATT > 0 is equivalent 
to 

A{e) (f (te) - f (oe)) ■ (f (fe) - f (oe)) > 0. 

e£E 

From now on, we assume, together with Q, that A{e) is symmetric 
and positive dehnite (this condition is natural if we assume that the 
interaction among atoms is a superposition of two-body interactions). 


Example (1) (Monoatomic crystal lattices) This is the case that Xq 
is a bouquet graph (i.e. the case that the unit cell contains exactly one 
atom). If A{e) is symmetric, then the condition © is satished since 
v(e) = —v{e). 

(2) (The scalar model) This is the case that A{e) = a{e)I with a 
positive-valued function a(e) on E such that a(e) = a(e). If 


(7) Y = 0’ 

e^Ex 

then the condition (jHI) is satished. It should be interesting to point out 
that ^ satisfying o is a discrete analogue of (vector-valued) harmonic 
functions, and that $ induces a “harmonic map” of Xq into the hat 
torus R^/L (see [T^L 


We conclude this section with a proof for the following lemma which 
we have employed to deduce the condition 
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Lemma 3.1. Let T G M„(R) be a real n x n matrix. Ifn>3, then 
there exist Ui,... ,Un G SO{n) and real scalars Ci,... ,cn such that 

T = c,{U,-I) + --- + cn{Un-I). 

Proof. It is enough to prove that, if tr T(t/ — J) = 0 for every 
U G SO{n), then T = O (the zero matrix). Take a skew-symmetric 
matrix S. Then e^^ G SO{n) for t G R, so that differentiating both 
sides of tr Te^^ = tr T, we hnd tr TS" = 0. Using 

tr T*S = tr ST* = -tr S*T* = -tr TS = 0, 

we obtain tr {T — T*)S = 0. Thus T = T*. Without loss of generality, 
we may assume T = diag(Ai,..., A„), a diagonal matrix. For U = 
{uij) G S'O(n), 

tr T = tr TU = Ai^n XnUnn, 


so that 


Ai A„ 

Applying this to 


( cos 6 
sin 6 

o 


— AiMii XnUnn- 

— sin 6 ^ O \ 
cos 9 , 

In-2 


we obtain 


Ai Xn — (Ai -|- A 2 ) cos 0 -|- A 3 Xn, 


which implies that Ai -|- A 2 = 0. In the same way, we have A, -|- Aj = 0 
for i 7 ^ j. Therefore we conclude that Ai = ■ ■ ■ = A„ = 0, and hence 
T = O as desired. □ 


Remark If n = 2, then the claim of Lemma o is not true. 


4. Continuum limit of a crystal lattice 

The condition © was deduced by a rather forcible and formal ar¬ 
gument. However, this condition turns out to be natural as we shall 
see below, once we agree that a crystal lattice is a discretization of a 
(uniform) elastic body, or the other way around, the continuum limit 
of a crystal lattice is an elastic body. We will also see that (P) plays a 
signihcant role in the later discussion. 

We shall infer what the continuum limit of the crystal lattice X 
should be by comparing lattice vibrations with elastic waves. The 
outcome is used to give a physical meaning to the constant cq, but is 
not required to establish the asymptotic of the density of states. 

We hrst recall that, in general, an elastic body is characterized by the 
mass density p and the elastic constant tensor Catyj satisfying Caipj = 
Ciafij = Caijfs = Cpjai (see [ini)- An elastic wave f = f(t, x), t G R, 
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X G R^, propagating in the elastic body, is a solution of the wave 
equation 



where, for each i, j, Aij is a matrix-valued function on dehned by 
= Caipj (the («,/3)-component of the matrix Aij). 

We easily see that ^Aij = Aji. 

We assume that the elastic body obtained as the continuum limit is 
uniform in the sense that p and Aij are constant. Symmetrizing Aij if 
necessary, we may assume Aij = Aji. 

It is natural, from the nature of mass density, to put p = m(Vo)/V. 
To surmise the form of the matrix Aij, take a smooth function f : 
R X R^ —> R^, and dehne : R x —> R^ by setting 

U{t,x) = 

(thus is a discretization of f with respect to the space variable). Then 


-{KU){t,x) 


A(e) [f(^t, S^{te)) - 5<^{oe))] 


di 




j=l e£Ex 
3 


1 02f 

9^' ^ 


i,j=i eeEx 


where v{e) = (n(e)i,n(e) 2 ,^( 6 ) 3 ). Under the condition 0. the hrst 
term vanishes. Therefore, we have, for a sequence {a;^} in V with 
lim^io = X, and a sequence in R with lim^jo = t, 


lim • 

(510 


-6 ‘^{Kfs){t5,xs) = -Y ^(e)*'^(e)i^(e) 


i,j=l eSEo,: 


dxidxj 


(^,x). 


provided that {a:^} is in the orbit containing x G JU. On the other 
hand, 

p2c pp.c 

\\m5-^m{x5)^{t5,X5) = m(x) —(t,x). 

Thus, taking the sum over the fundamental set JF, and dividing by V, 
we may presume that the equation of motion approaches the equation 


dH 


1 


3 

i,j=i eei?o 


dH 

dxidxj 
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as the mesh of the lattice becomes hner. This implies that the (sym¬ 
metrized) elastic constant tensor of the elastic body corresponding to 
our crystal lattice is given by 

^ v{e)iv{e)jA{e). 

eS-Eo 


Remark The differential operator (the elastic Laplacian) 



d 

dxj 


can be written as —d*d, where 

d : A°(R3, R3) —. yl^(R3, R3) 

is the exterior differentiation acting on R^-valued differential forms, 
that is, 

3 

i=l 

The operator d* is the formal adjoint of d with respect to the inner 
products on A*(R^, R^) {i = 0,1) dehned by 


{i,g)= f f'gprfx (f,geT(RlR=>)), 

v)= ■ Vi V e A^(R^, R^)), 

where u = (cti, ct 2, c^s), rj = (?7i, 172, i/s)- 


We now return to the difference operator D. We easily see that D 
is an L-equivariant linear operator of C(V, C^), the space of C^-valued 
functions on V. If we dehne the Hilbert space i‘^{V,m) by 

m) = {f G C{V, C^); ||f |p ;= 52 ‘ ^{x)m{x) < cx)}, 

x£V 

then D restricted to the subspace i'^(y,m) is a bounded self-adjoint 
operator of i‘^(y,m). We call D the discrete elastic Laplacian. This 
naming is justihed by the expression 


D = -d*d, 

where d : i'^(y,m) i‘^{E,A), a discrete analogue of the exterior 

differentiation (or the coboundary operator in cohomology theory), is 
dehned by 


di{e) = f(te) — f(oe). 
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and 


\E,A) = , 7 (e) = -r/(e), 

]-^A{e)7]{e) ■ 7]{e) < 00 }. 


2 


e£E 


Indeed, the explicit expression for d* is given by 

1 


(9) 


{d*u}){x) = 


m{x) 


^(eV(e)- 


^■GEx 


This is checked by the following computation. 

{di,u;) = ^^v4(e)(f(te)-f(oe)) ■o;(e) 


eeE 

e&E 

^ f(^e) ■ A{e)uj{e) - ^ ^i{oe) ■ A{e)uj{e) 

eGE 


eGE 


^f(oe) ■kl(eV(e), 


eGE 


where we have used the assumption ^A{e) = A{e). The last term is 
written as 

1 


Ef( 

xGV 


X) 


m{x) 


kl(e)o;(e) jm(x), 




from which © follows. 

Therefore the operator D is not only a discretization of the elastic 
Laplacian, but also its conceptual analogue. 


5. Hamiltonian formalism for lattice vibrations 

This section is devoted to a brief explanation for the quantization of 
lattice vibrations which is performed in a slightly different way from 
the current one in the physical literature. To this end, we shall start 
with the Hamiltonian formalism for lattice vibrations. 

We put S = and denote by (•, ■) the inner product on the 

Hilbert space S. We shall regard S' as a symplectic vector space with 
the symplectic form ut dehned by 

uj{u,v) = lm{u,v) {u,v ^ S) 

(Im 2; denotes the imaginary part of z G C). Dehne the Hamiltonian 
Hby 

H{u) = ^{V-Du,u). 
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Then the Hamiltonian equation for H is given by 


du 

dt 


—^/—iV—Du, 


which is obviously equivalent to the equation 


dH 

dt‘^ 


Di. (The equiva¬ 


lence can be established, for example, by the relation f = Ren). 

We wish to quantize the Hamiltonian system (S', a;, H). To avoid the 
difficulty arising from the infinite-dimensionality of S, we shall decom¬ 
pose {S,u},H) into a direct integral of finite dimensional Hamiltonian 
systems. The idea, which essentially dates back to Bloch’s work on 
periodic Schrodinger operators and has been taken up in a different 
manner by physicists, is to use the irreducible decomposition of the 
regular representation of L on the Hilbert space 

e{L) = {f:L^C- 


( 10 ) 



(x,C)dx, 


where L is the unitary character group of L, and dx denotes the Haar 
measure on L, which is normalized so that 


Jdx=l. 

Similarly to m, we may construct a direct integral decomposition 


/*© 

(11) {S,u;,H) = {S^,u;^,H^)dx. 

Namely, we take 


S^ = {u:V ^ C^; u{ax) = x{(^)u{x)}, 
uj^{u,v) = Im(M,n)x, 


^x(“) = 


where 

(M,n)x = E u{x) ■ v{x)m{x) 

x^T 

is the scalar product in 5^. Note that dim Sy^ = 3n, and Sy is an 
invariant subspace (in fact, even an eigenspace) of the action of L in 
C(H, C^) (the space of all C^-valued fuctions on V\ Late we will 
sometimes refer to the elements of Sy as Bloch functions, and to the 
components of different objects in the direct integral decomposition 
above as the Bloch components. 

Note that the action of L in Sy is not irreducible (unlike the corre¬ 
sponding subspace for the action of L in The isometry between 
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S and the direct integral 

r® 

is given as follows. For u E S with hnite snpport, dehne G by 

( 12 ) u^{x) = 

a^L 

Then the extension of the correspondence u i—*• {mxI gives rise to the 
desired isometry. The inverse isometry is 


(13) e X e t| ^ m(x) = 

where we assnme that the fnnction % i—^ is sqnare-integrable on L. 

The “twisted” operator is dehned to be the restriction of D : 

C{V,C^) —^ C{V,C^) to C C{V,C^). It should be pointed out 
that, in physical terms, the direct integral decomposition (HH) corre¬ 
sponds to the “sum” over wavenumber vectors. 

Since d and d* are also decomposed as 


r® 

d = dx dx, 

d” = r(d’)x dx, 

and (d*)^ = {d^)*, we have and hence < 0. Fur¬ 

thermore, if X 7 ^ 1 (1 being the trivial character), then < 0. Note 
that 0 is an eigenvalue of —D^ of multiplicity three whose eigenfunc¬ 
tions are constant, because we assumed X to be connected (hence Xq 
is connected too). We enumerate the eigenvalues of —D^ as 

0 < Ai(x) < A2 (x) < • ■ ■ < A3„(x). 

The functions Aj(x) are continuous on L, and the hrst three eigenvalues 
Ai(x)) A 2 (x )5 A 3 (x) are perturbations of the eigenvalue 0 = Ai(l) = 
A 2 ( 1 ) = A 3 ( 1 ) of —Di, which are said to be acoustic branches, while 
other eigenvalues are said to be the optical branches. 

Remark. The spectrum a{—D) coincides with 

3n 

IJ Image Afc. 

k=l 

Therefore —D has a band spectrum in the sense that a{—D) is a union 
of hnitely many closed intervals. In the case of monoatomic lattices, 
we only have acoustic branches, and a{—D) is an interval. 

Choose an orthonormal basis ei,..., 63 ^ of such that 

Aj(x)ej. 
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Then the R-basis {ei,..., e^n, , y/^esn} yields a canonical 

linear coordinate system (pi,... ,p 3 n, Qi, ■ ■ ■, Qsn) by setting 


3n 

U = '^{qi + y/^Pijei. 
i=l 

In terms of this coordinate system, the Hamiltonian is expressed as 

^ 3n 

= 2 ^VHx){p^i + q^)- 

i=l 


Note that = 2 + ^i) is the Hamiltonian of the sim¬ 

ple harmonic oscillator with the frequency A/Aj(x)/27r. Therefore the 
quantized Hamiltonian for is given by 


1 3n ^_ 02 3n ^ 

= 2 E yvb( - + 

i=l i=l 

The standard fact tells that the spectrum of the quantized Hamiltonian 
^x,i (t 7^ 1) consists of simple eigenvalues 

= WMx) (k + (A; = 0,1,2,...). 

The “quantum particle” with the energy -ExlA be a phonon. 


6. Computation of the internal energy 


Following a general recipe in quantum statistical mechanics (more 
precisely, the Bose-Einstein statistics for phonons), we dehne the par¬ 
tition function for H^ i by 


k=0 


^-fix/x!X^/2KT 
1 _ ^-tiy/X^IKT' 


Then the internal energy of for the Gibbs distribution 
Pk = {k = 0,1, 2,...) 

is given by 


UxAT) = ^ E^,i,kPk = hx/\i{x) 


k=0 


1 1 
2 ^hx/xX^)IKT _ ^ 


In view of the additive property of the internal energy, it is natural 
to dehne the internal energy of the lattice vibration (per unit cell) by 
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which is expressed as 

I /, tr dx + /. tr A = C/„ + C/.(T). 

To transform this formula into one in terms of the density of states, 
we employ the notion of von Neumann trace (or L-trace; see PP), which 
is defined, for a L-equivariant bounded linear operator T : S' —> S', as 

tr£, T = ^^tr t{x,x)m{x), 

where t{x,y) is the kernel function (matrix) of T, that is, 

(14) T{{x) = E t{x,y)^{y)m{y). 

y&V 

Note that t{x, i/) is a 3 x 3 complex matrix for fixed x and y, and the 
matrix-valued function x,y ^ t{x, y) is invariant under the diagonal 
action oi L onV xV given by o'{x,y) = {ax, ay). In particular, the 
scalar function x i—> tr t{x, x)m{x) is L-invariant on V, hence the trace 
tr^ T in the definition above does not depend upon the choice of the 
fundamental domain T. 

We should point out that the notion of L-trace is introduced in a 
more general setting, say in von Neumann algebras of type IIi. 

The L-trace is related to the direct integral decomposition. The fol¬ 
lowing Lemma is known in the theory of operator algebras (in a more 
general context), but we will supply an elementary proof for conve¬ 
nience of the reader. 

Lemma 6.1. (Trace Decomposition Formula) Let T : S ^ S be a 
bounded linear operator commuting with the action of L, and 

r® 

is its direct integral decomposition. Then 

(15) tr^ T = y tr dy. 

Proof. First we will introduce some notations. Let us identify the uni¬ 
tary character group L with the torus 

Jl = rVl* 

via the correspondence 

y e ^ y e L, 

where 

X{a) = exp (271 ^^y ■ a), 

so that a small neighborhood U{1) of 1 in L is identified with a neigh¬ 
borhood t/(0) of 0 in R^. 
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We will present T in terms of operators 

T°:C'(l/o,C=')^C'(W,C='), 

where CiV^, C^) is the Hilbert space of all C^-valued functions on Vq, 
with the scalar product 

(n,n) = E u{x) ■ v{x)m{x), 

x£Vo 

dimcC(Vo,C^) = 3n. The operator T° is defined as T° = 
where : C(Vo, C^) —>• is a unitary operator, given by 

(UX){x) = s^(x) = 

where x E V, 7i{x) is the canonical projection oi x to Vq = V/L. It is 
easy to check that G S^. Indeed, 

s-^{ax) = exp (27rv^^x ■ ^(c’'2;))f(7r((ja;)) 

= exp ( 27 rA/^x ■ (^(3^) + c’'))f(^(2^)) 

= exp{27r\/^x- 
as desired. The inverse operator 

y-‘: ^ C{V„, C’) 

is given by 

where x E V. (The left hand side is well dehned because the right hand 
side is easily seen to be L-invariant.) 

The advantage of using the operators T° (compared with T^) is that 

they act in the same space C{Vo, C^) for all y G L. 

Now let us present the operator T° by its kernel function (matrix) 
= t^{x, y), where x,y E Vq, so that 

(T°f)(x) = ^fJ(T,j/)f(j/)m(j/). 

y&Vo 

Let us choose a function u E S with a hnite support, and hnd its 
Bloch components (y G L), dehned by (IT^ . which can be also 
rewritten in the form 


u^{x) = ^’^'^^■°'n(cra;). 

(jGL 


The Bloch components of Tu are {Tu)^ = ^ u^. Let us 

calculate them explicitly. We have 

{U~^u^{x) = u^{x) = 

aSL 


E 


= > e 




u{ax), 


(jGL 
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which is L-invariant (hence well defined on Vq). Further, we find 

y(iT,a&L 

y&V 

Therefore, 

{T^u^){x) = (f/,T°f/,-^J(vr(x)) 

= X] ^5(7r(a;), 

y&V 

Taking into account (HI, we see that 
{Tu){x) = jjT^u^){x)dx 

This means that the kernel function of T is given by the formula 

( 16 ) t(x,y) = /e=-'^’<'«*>-*«>i;(i(i),,r(j,))<ix. 

Note that by the general properties of direct integrals 
||r : S' ^ S\\ = ess sup{||Tx||; y e L} = ess sup{||T°||; y G L}. 

In particular, T is bounded if and only if the norms ||T^|| = ||T°|| are 

a.e. uniformly bounded for y G L, or, equivalently, the matrix elements 
(x, y) are uniformly bounded a.e. in y for all x,y ^ Vq. 

Now taking y = x in dH, applying matrix trace (to the 3x3 matrices 
in both sides) and summing over x G JF, we obtain 

tr^ T = ^tr t(x,x) = / ^ tr (7r(x),7r(x))(iy 
xe.u 

= [ tr T^dx = [ ti T^dx, 

Jl Jl 

which ends the proof. □ 

Remark. The kernel representation formula (unD also allows to provide 
a relation between the smoothness of the function y i—>■ T° and the off- 

X. 

diagonal decay the kernel t{x,y), i.e. its decay as \${x) — ^{y)\ —>■ cxo. 
More precisely, let us say that a bounded linear operator T : S' —> S 
(commuting with the action of L) is smooth, if the matrix elements 
t^{x,y) are in C°°{L) with respect to y, for any x,y E Vq. If T is 
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smooth, then, integrating by parts, we can for any positive integer N 
rewrite (CHD as follows: 

t{x,y) = (1 + \^{x) X 

j^^2nV^x-mx)-<P(y)) [("x _ ^Ax)^f°(7r(a:),7r(i/))] dx, 
where denotes the Laplacian with respect to y. It follows, that 
(17) \t{x, y) I < ^^(1 + \<d>{x) - <d>{y) A = 0,1, 2,.... 

Vice versa, if these estimates hold, then we can calculate T° as the 
operator U~^TU^, identifying elements of (7(1/0, C^) with L-invariant 
C^-valued functions f on V, and then applying T to the Bloch function 
f/^f using the kernel (matrix) representation of T in terms of t{x,y) 
(see (HI). This leads to the expression 


{U-^TUJ){x) = i/)m(i/)f(i/) 

y&V 


_(tGL 


m{y)^{y), 


where the sums converge due to m- The latter expression is equiv¬ 
alent to the following presentation (which is the inversion formula to 
(fT^) for the matrix : 


(18) tl{Ti{x),Ti{y)) = '^e ‘^(^y'))t{x, ay), 

(jGL 


where the left hand side is well dehned because the right hand side is 
L X L-invariant on V x V, i.e L-invariant separately in x and y. We can 
also differentiate (HHD with respect to y termwise , which implies the 
smoothness of the function y i—>■ T°. □ 


Remark. From m we also deduce 

tr = tr T° = ^ ^ g-2,rFTTx-(<^(A-<^(-A)tr t{x, ax), 

x^T (tGL 

which is an analogue of the trace formula (see uni). 


Let 


-D 


A dE{\) 


be the spectral resolution of —D. Then the integrated density of states 
is dehned by 

(p{X) = tiL E{X). 


We hnd that (p is a nondecreasing function, and 
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SO that 


dip{X) = 3n. 


We should note that, when we take a decreasing sequence of lattices 


L = Lo D hi D La D ■ ■ ■ 

with n^ofunction (^(A) coincides with the limit 

lim [#(L/Li)]“V*(A) 

1^00 

at continuity points of (p, where (pi is the counting function of eigen¬ 
values for the periodic boundary value problem: 

-Df = Af, 

i{ax) = i{x) {a E Li) 

(see [inj, [HI), that is, if Ai < ■■■ < {N = ^{V/Li)) are the 
eigenvalues for the boundary problem above, then V3j(A) = #{Aj < A}. 
Thus the function (p deserves to be called the (integrated) density of 
states. 

Now using 


(19) 

we get 


f(-D)= JJ(-D^) dx, 
j ti: f(-Dy,)dx = tri /(--D) = j 




U,(T) = 


dip{X), 

/■“ hVx 

Iq ^hVX/KT _ 1 


dV3(A), 


C{T) = 


dUi 

W 


roc f?\ Hx/X/KT 

/ _ 

/o [e^VXlKT _ 1) 


jd(p(A). 


Uo is what we call the zero-point energy^ which indicates that vibrations 
in a solid persist even at the absolute zero of temperature, a peculiar 
feature of quantum mechanics. 


Remark. (1) It is straightforward to see that 

supp dip = spectrum(—D). 

(2) We may consider integrated densities of states for other bound¬ 
ary conditions. From the abelian nature of the crystal lattice, however, 
it follows that the notion of density of states is independent of the 
particular boundary conditions imposed. 
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(3) The partition function for the Hamiltonian H^ i in the classical 
statistical mechanics is given by 


ZxAT) = 


exp - 


/ exp 

'R2 

2Ktt 


KT 

2KT 


dpidqi 


(p'i+Qi)) dpidq, 


--T, 


and hence the internal energy is computed as 


U^AT) = A'T. 

From this, it follows that the internal energy of the solid (per the unit 
cell) in the classical setting is given by 

„ 3n 

LJ2^xAT) dx = 3nKT, 

i=l 

which is nothing but the law of Dulong-Petit. 


7. Asymptotics of Lp {\) at A = 0 

Our task is now to establish asymptotics of Ui{T) and C{T) as T 
goes to zero. To this end, we need to study the asymptotics of 93 (A) as 
A i 0. 

If we put 

= #{b^i(x) < A}, 

then, due to (USD, 

pW = dx, 

which, as is easily checked, is equal to 

vol({x e L; (px(A) = 1}) + 2 vol({x € L; (^^(A) = 2 }) 

+3vol({x e L; (p^(A) = 3}) 

= vol({x e L; Ai(x) < A < A 2 (x)}) 

+2vol({x e L; \ 2 {x) < A < A3(x)}) 

+3vol({x e L; A3 (x) < A}) 

3 

= ^vol({x; A„(x) < A}) 

for sufficiently small A > 0. Hence it suffices to establish the asymp¬ 
totics of 


vol({x; A„(x) < A}) (a = 1,2, 3). 
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We will use the notations introduced in the previous section and start 
with the operator D. Note that D has a “finite radius of interaction”, 
that is, its kernel function d{x,y) vanishes if \^{x) —^{y)\ > R, where 
i? > 0 is a constant (which can be taken to be maxe eE{\^{te)-<P{oe)\}). 
In particular, D can be naturally extended to all C^-valued functions 
on V. This allows to decompose D into direct integral of the Bloch 
components explicitly by taking to be the restriction of D on 

Let us also explicitly calculate the operator 

Dl:CiVo,C^)^CiVo,C^). 

We claim that 


Indeed, let us first move to the Bloch space via the correspondence 
f e C{Vo, C^) 1-^ s(x) = exp (27r\/^x ■ lP(a;))f(7r(a;)), 
and take 

S' 


_g27rv^X-‘^(oe) 


f(7r(oe))] 




X—- y - f(7r(oe))l 

mix' 


eeE^ 


— g27rv^X-‘^ W . 


m(x} 


y kl(e) (te) - f (oe)]. 


e^Ei 


071(2:) 


Dividing by we get back to ^(Vo, C^) and obtain the de¬ 

sired expression for D^. 

Let us take y G R^. Since is a one-parameter family of symmet¬ 
ric matrices depending analytically on t G R (where |t| is small enough 
so that tx G t/(0)), one can apply the following result (see jS]). 


Theorem 7.1. Let A{t) be an analytic family of symmetric matrices 
of degree n. Then there exist analytic functions Xi(t ),..., Xn(t) {called 
analytic branches of eigenvalues) and orthonormal basis ui{t),, Un{t) 
depending analytically on t such that 

A{t)ui{f) = Xi{t)ui{t). 


The acoustic branches Aq(x) are not necessarily analytic in x because 
they are possibly bifurcated, not only at x = 1, but also at other x- The 
theorem above, however, says that one can find e > 0 such that Xa{tx) 
{a = 1, 2, 3) are real analytic in t G [0, e). To see this, take continuous 
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branches , ^ 2 (t), fisit) of eigenvalues of with = 0, and 

let lij = {t G R; /ij(f) = {i ^ j). Then I^j is either discrete or 

R. Put 

e = niin{|f|; t ^ 0, t E hj j}. 

If 0 < t < e, then ^ or /Xj = /i^, so that are an¬ 

alytic on [0, e) for all j = 1,2,3, and Xa(tx) coincides with one of 
iJ, 2 {t), Hsit) on [0,e). 

Since > 0 for all t G R, j = 1, 2, 3, we conclude that 

X 

SaixY ■= lip T^Xaitx) (Saix) > 0) 

HO 

exists. The quantity 

(27r||x||)-^s„(x) 

is said to be the acoustic phase velocity for the direction y. It turns out 
to be the phase velocity of elastic waves in the uniform elastic body 
corresponding to the crystal lattice (see Section HDD. 

Recall that we assume Ai(x) < A 2 (x) < A 3 (x). From the dehnition 
of Sa(x) we deduce that si(x) < S 2 {x) ^ 'S 3 (x)) ^ind it follows from 
Theorem o below that Sq,(x) is continuous in y. We also have 

Sa{tx) = tsaix) {t>0). 

A perturbation argument leads to 

Theorem 7.2. 1) Sq,(x)^ {a = 1,2,3) are eigenvalues of the symmetric 
matrix 

^ ' e£Eo 

In particular, Sa(x) > 0 for x 7 ^ 0, and 

27r2 

Si(x)^ + S 2 {x? + S 3 (x)^ = tr ix ■ v{e)fti A(e). 

2) Sa(x) are piecewise analytic m y G R^ in the sense that there 
exists a proper real-analytic subset S in R^ such that s^ is analytic on 

R\S. 

3) The Lebesgue measure of the set {x G R^; So(x) = 1} is zero. 

Proof We put Dt = and Xa{t) = Xa{tx). In view of Theo¬ 
rem o one can hnd an orthonormal system fi,t, f 2 , 4 , fs.t £ C{Vo,C^) 
depending analytically on f > 0 for small t, so that 

(20) -Afa,i = A„(t)f„,i (a = 1,2,3). 

Recall that is a constant function. If we dehne ea{x) by 

e«(x) = 

then we observe that {ei(x), e 2 (x), e 3 (x)} forms an orthonormal basis 
of CT 
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Differentiating both sides of (EHD with respect to t at t = 0, we obtain 

-^ X] A{e) [2'Ky/^{x ■ u(e))fa,o(ie) + fa,o(te) - f„,o(oe)] 

eSEo. 

= A„(0)f„,o(2:) + A„(0)fo,o(2:)- 

Since fc^o is constant, and Aq(0) = Aq( 0) = 0, by using Q, we have 

'Y' A(e)(fa,o(te) - fa,o(oe)) = 0 

mix) 

so that fc^o is constant. 

Differentiating (EOl) twice, we obtain 

- ^ Y Al(e)[-47r2(x-n(e))2f„,o(te) 

mix) 

+471 • n(e)f«,o(te) + fa,o(^e) - fa,o(oe)] 

= A„(0)f«,o(a:)- 

Using the fact that is constant, and again ©, we have 

-^ Y [ ~ ■ t'(e))^f„,o(te) + fa,o(te) - f„,o(oe)] 

7Y}/ \X) 

= A„(0)f„,o(2:)- 

Here taking the inner product with f^ o; and noting 
(-Dfa,05 f/3,o) = (fo,0 5-Df/3,o) = 0, 

we hnd 

4vr^ Y (a ■ v{e)YA{e)ia^o{te) ■ f/3,o(oe) 
eSEo 

= A„(0) Yj ^a,o{x) ■ if)^o{x)ni{x) 
xsVb 

Aq,(0)(5q.^ 2Sq.(x) Saf3i 

or, equivalently, 

Alxea(x) • epix) = SaixfSap, 

and hence Sa{xY i® an eigenvalue of with the eigenvector eQ(x). 
Since we assumed that A{e) is strictly positive for every e, we obtain 
also that Sq(x) > 0, a = 1, 2, 3. 

To show 2), let D = D(x) be the discriminant of the cubic polyno¬ 
mial 

det(z/ - A^) = + b{x)z^ + c{x)z + d{x)- 

If D ^ 0, then Sq,(x) (a = 1, 2, 3) are distinct and analytic on R^\5o 
where 


5o = {x e D(x) = 0}. 
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In the case D = 0, define the polynomial Di(x) by 
Di(x) = ‘^Kxf - 9&(x)c(x) + 27d(x). 

Note that the equation det( 2 :/ — A^) = 0 has a root of multiplicity one 
(hence another root has multiplicity two) if and only if Di(a) ^ 0. If 
Di = 0 , then 

= S2(x)^ = ss(x)^ = -h(x) = 

hence is analytic on R^\{0}. In the case Di ^ 0, define a proper 
analytic subset iSi by 

5i = {xGR^ Di(x) = 0}. 

One can easily check that 6 (x)^ — 3c(x) 7 ^ 0 on R^\{0} and the roots 
of det( 2 ;/ — Ay^) = 0 are 

9d(x) - &(x)c(x) 


2(6(x)2-3c(x)) 

-Kxf + 4fe(x)c(x) - 9d(x) 


(multiplicity two), 
(multiplicity one), 


Kxf - 3c(x) 
so that Sq is analytic on R^\5i. 

The assertion 3) easily follows from 2), if we take into account the 
explicit form of and strict positivity of A{e). 

This completes the proof of the theorem. □ 

Think of Jl as the fiat torus with the fiat metric induced from the 
Euclidean metric on R^. The normalized Haar measure on L is identi¬ 
fied with 

vol(JL)“^dx, 

where dy, in turn, denotes the Lebesgue measure on the Euclidean 
space R^. Note that V = vol(J 2 ,)“^ is the volume of a unit cell. 

Put 

^a(A) = {x G R^; Aq,(x) < A}, 

and denote by the characteristic function for a subset A in R^. Then 


vol({x e L- A„(x) < A}) = V 


'Ac(A) 


dx = ^ 1 a ,( a )( x ) dx 


'R3 


VA’/" / 1.4.(A)(^\') dx. 


'R3 


Since 


1ama)(^x) = 


1 if A ^XaiVXx) < 1 
0 otherwise, 

we observe, due to Theorem iia that 


hmlA4A)(^x) = 1 a,(x) 

A].U 


a.e. 
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where 

Aa = {xe Sa(x) < !}• 

From this, it follows that 

hmvol({x e L] A«(x) < 

A].0 

= V / Ujx) dx- 

Jr 3 

Using the polar coordinates (r, hi) G R+ x we obtain 


3 Js2 


[ Iac{x) dx= [ r'^drdVL = \ [ — dQ. 

d{{r,n); Sc(H)r<l} 

In conclusion, we have 

Theorem 7.3. 


3 

= iv / V 




dQ. 


Remark. The asymptotic formula for an acoustic branch 

VAo(rx) ~ rsaix) (r i 0) 
is said to be the linear dispersion law. 

8. Rigorous derivation of the T^-law 
We are now ready to prove the T^-law. Put 


( 21 ) 




3 Js^ Sa{^y 
We have shown in the previous section 

99 (A) ~ CqA^/^ as A I 0. 

Making the change of variables x = hy/X/KT in the integral 

nVx 


Ui(T) = 


^hV\/KT _ 


d^{\), 


we obtain 


/■oo ^ 

C / i ( T ) = fi - VCT ' / .^—(wt)-^Ap((wtxX), 

Jo e — i 

where wt = KT/h. Using 99 (A) we have 

\im{wT)~^(p{{wTxy^ = cqx^, 
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and hence 


r-oo ^3 


-1 


lim[/i(T)T-^ = 3H-^K^Co 
Tio 

5 

Here we have used the well-known identity 

TT^ TT*^ 


dx 


Similarly we have 

C{T) = fr^K^T^ 

hence 


r-oo 


(e* - 1)2 

lirncmr-^ = 3r^ir% f 

no 7 o 


(wt) ^d(p{{wTxf), 


oc 


0 (e"-l) 


dx = -ti'^coH ^K^Cq 

Z R 


since 

:- ^dx = 4!C(4) 

gX _ 1)2 '’V y 

Summarizing our computation, we have 

Theorem 8.1. Hs T J, 0, 




17i(T) ~ 

5 

cm ~ -Tr'^coH-mm. 

5 


9. Isotropic case 


In order to compare our result with Debye’s asymptotic formula, we 
suppose that the continuum limit of the crystal lattice is isotropic, that 
is, the symmetrized elastic constant tensor 

eeEo 


satishes 

3 

'^{Aj)apXiXj = (a + b)XaX0 + bdafsWxf, 
ij=l 

where y = (xi,X2,X3)5 and a,b are positive constants which are what 
we call Lame’s constants in the theory of elastic bodies. Recalling 
which was introduced in Theorem o we obtain 


3 

A = 4Ap~\iXjAj {p = m(Ro)/V). 
*j=i 
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The eigenvalues of are 

47r^(a + 26)p“^||x||^ (multiplicity one), 
A7r^bp~^\\x\\^ (multiplicity two), 


since 


and if X ■ X 


Therefore 


3 

X](^x)«/3X/3 = 471^(0 + 2b)p-^\\x\\‘^Xa, 

/9=1 

0 for X = {xi, X2, X 3 ), then 
3 

^{A^)c,0Xp = ATT^p-^WxfXa. 

/3=1 


( 22 ) 


si(x) 


S2(x) = 271 


y^llxll, Ssix) 



' a + 26 . 


According to the theory of elastic waves, we say that ci = \ i -is 

the longitudinal phase velocity, and Ct = \ - is the transverse phase 

Vp 

velocity. In terms of these phase velocities, we have, by substituting 
(1221) for (EH) 


V 


Co = 


(— — 

Gtt^ Vc? c? 


We thus recover Debye’s result for the isotropic case. 


10. Final remarks 


In general, the “plane wave” 

f(t, x) = exp [v^(27rx ■ y ± tSa(x))] ea(x) 
in a uniform elastic body is a solution of the elastic wave equation 

^ _ A d^f 

^ dt"^ ^ dxidxj' 

i,j=i ■> 


To get real solutions we can take real or imaginary parts (or their real 
linear combinations). For example, 

f(t, x) = cos (27rx ■ X ± ts„(x))ea(x) 

is a real plane wave which solves the same equation. 

It is interesting to point out that the “discrete plane wave” 


i^{t,x) = exp 


\/^(2vrx ■ ^{x) ± tV'Ac,(x)) 


gx(^(^)) 
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in the crystal lattice is a solution of the equation 


dH 

dt‘^ 


= Di 


provided that = AQ,(x)g^. (The real discrete plane waves can be 

again obtained by taking linear combinations of the real and imaginary 
parts, which do not, however, have the cosine or sine forms because the 
operator and the eigenvectors g^ are not generally real.) 

We easily observe that, for an acoustic branch Aq,(x) with lim^j^o g<5x = 
eQ,(x), we have 


= exp [v^(27rx-x±ts„(x))] e„(x). 


where the sequences and are supposed to satisfy hm^^o = t 
and lim^jo d<d>{xs) = x. This justihes the statement that the lattice 
vibrations approch elastic waves in the continuum limit. 

Finally, we shall check, using the von Neumann trace again, that 
the function coincides with the integrated density of states for 

the elastic waves (thus justifying Debye’s observation in his continuum 
theory). For this, we put 


V 





dxidxj ’ 


and let K{t, x, y) be the kernel function of the operator namely it 
is the fundamental solution of the parabolic equation 


9f 

'Wt 


= Vf. 


Then the L-trace of e® is given as 


tr^q 


w 


= j tr iF(t, X, x) dx. 


where P is a unit cell. We can easily show, by using the Fourier trans¬ 
formation. 


P(f,x,y) = (27r)-3 [ dy. 


3 

Mx) = p~^ AiXiAj 




where 
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Therefore 

tr^e® 

n i*00 

= V tr i^(t,x,x) = (27r)-=^V dil tr dr 

Js^ Jo 

n n QQ 3 4. 2 

= {2tt)~^V J ^dil J r^^exp dr. 

Using the equality 

^2g-ar2 ^ !1 ^q-3/2 o > 0, 

4 

we get 

tr^e® = /■ \2sa{^)~^ dn. 

If we denote by 9?o(-^) fhe integrated density of states for elastic waves, 
that is, L-trace of the projections in the spectral resolution of V, then 

d(foW = trre*^ = ^y/7rcot~^^^, 

so that, taking the inverse Laplace transform, we obtain 

3 

(^o(A) = dn = coA^/^ 

POO 2 

as desired (recall the identity / d(A^/^) = -y/7rt~^^‘^). 

Jo 4 
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